HZ''''£1'‘SYHNW2S8I 


2 Af. Encke , Orc a 2Vew> Method of 

MASSILIA. 

Elements . By M. G. Rumker. 
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These elements represent the middle observation, thus : 

Computation — Observation j +°!+ “ 

An ephemeris based on these elements was sent by M. G. Rumker at the 
same time: this has be.en forwarded to our principal observers. 


The observations of the solar eclipse of July 28, 1851, which 
form Part I. of Vol. XXI. of the Memoirs , are now published. 
There are six plates and 120 pages of letterpress. The price to 
Fellows is 5 s., and 105. to the public. There is a separate title- 
page for the purchasers of the part who may not want the volume. 


On a New Method of Computing the Perturbations of Planets . 

By M. Encke, Director of the Royal Observatory of Berlin. 

Translated into English and Illustrated with Notes by G. B. 

Airy, Esq., Astronomer Royal. (Appendix to the Nautical 

Almanac for 1856.) 

The main object of the author is to supply an easy and accu¬ 
rate method of computing the perturbations of the various small 
planets revolving between Mars and Jupiter , the number of which 
is constantly increasing from fresh discoveries. The method, how¬ 
ever, is applicable to any of the other bodies of the planetary 
system. 

Let t denote the time from which it is proposed to compute the 
action of the disturbing body. By a slight transformation the dif¬ 
ferential equations of the motion of the disturbed planet are reduced 
to the following form :— 

*=ff( Y + £(3 r- «))«*<*- 

?= #( Z+ 5< 
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Computing the Perturbations of Planets . 3 

#°, I/ 0 , 2 0 , being the elliptical co-ordinates of the disturbed planet, 
and r° its radius vector, as deducible from the elements of the 
orbit in which it was revolving at the instant t ; X, Y, Z, the re¬ 
solved values of the disturbing force in the directions of corre¬ 
sponding rectangular axes; jj, £, the perturbations of x , y , 2, 
counted from £; $r, the perturbation of the radius vector; and £ 
a constant quantity. 

These equations are integrated by mechanical quadratures, the 
successive numerical values of the quantity under the sign of the 
double integral in each equation being corrected by a peculiar 
process, founded on the theory of finite differences. 

The first step is to divide the period for which the perturbations 
are to be computed, into a number of sufficiently small intervals of 
time, each equal to r, and to calculate the values of X, Y, Z, for the 
times t — \ r, t + \ r, t -Jr, &c., using for this purpose the ellip¬ 
tical co-ordinates x°, y° , 2 0 , of the disturbed planet. As g, 

are assumed severally equal to zero at the commence¬ 
ment of the period for which the perturbations are to be computed, 
the first three values of the quantities under the sign of the double 
d % n d 2 Z 

integral (or — 2> ai ’ e readily found with sufficient accuracy 

by an approximative process, and from them the corresponding values 
of u, £, are deduced. For all the subsequent values of |, »),£, the 
computation is uniform. The mode of procedure is this :—From 

the three last values of ~f an approximate value of the same 

quantity for the following argument is readily deduced, and a pro¬ 
visional value of | is then found by means of the formula of inte- 
• d 2 £ 

gration | = "f-\- "f being a numerical quantity computed 

with the utmost facility by the aid of results already obtained. 
The corresponding values of n and £ being found by a similar 
process, the three resulting numbers are employed in computing 
the second part of the expression under the integral sign in 
each of the fundamental equations. Finally, the numbers thus 
found, when applied to X, Y, Z, give the corrected values of 
d 2 £ d 2 y d 2 Z . . . : . 

d¥ y dU y dl 2 ’ m ^eh' turn lead by double summation to the 

true values of g, »j, £. By a similar process the values of g, n, £, for 
the next argument are ascertained, and so on in succession, the 
operation being rapidly carried on throughout the entire period for 
which the perturbations are to be computed. 

In order to test the practical utility of this method, M. Encke 
compared its results with those obtained by the method of the 
.variation of arbitrary constants. For this purpose the perturba- 
tions of Vesta by Jupiter from 1854, May 24, to 1855, October 7, 
were subjected to a double computation. The following table 
exhibits the values of the polar co-ordinates of the disturbed planet 
at the end of this period, as resulting from the use of both 
methods :— 
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M. Encke, On Computing the Perturbations of Planets . 


Variation of Constants 

O / B 

Longitude 335 20 6’4 

Latitude — 5 40 42*3 

Log. Rad. Vector . 0*3665032 

M. Encke mentions as a proof of the easy application of this 
method, that in computing the perturbations of Vesta with given 
places of Jupiter for thirteen epochs, only eleven hours were 
occupied, notwithstanding that he was frequently interrupted in 
the course of the operation. 

If extreme accuracy is required, the values of x, y, z , resulting 
from the application of g, tj, £, to x° , y° , z°, may be employed in 
computing X, Y, Z, and then by repeating the operation we may 
obtain corrected values of g, »?, This process, which gives the 
perturbations depending on the square of the disturbing force, 
may be repeated over and over again until the assumed and re¬ 
sulting values of jj, £, finally agree. The method, however, is so 
rapidly exhaustive that this object is generally secured at the first 
repetition. 

When the perturbations are large the method is still practi¬ 
cable, but it becomes necessary in this case to reduce in a corre¬ 
sponding degree the magnitude of the intervals of quadrature. 

The following are some of the principal advantages attending 
this method of computation : — 

1. It proceeds from the fundamental equations of motion by 
a series of arithmetical operations of extreme simplicity, only a few 
analytical expressions of by no means a complicated character 
being employed in the computation. 

2. It affords the means of immediately correcting the quan¬ 
tities employed in computing the disturbing forces for each suc¬ 
cessive interval, and thereby leads rapidly to results of great 
accuracy. 

3. Since the perturbations of x, y,z , are of the same order as 
those of the polar co-ordinates, they serve to afford unequivocal 
indication of the extent to which the intervals of quadrature ought 
to be diminished. It is well known that the method of the varia¬ 
tion of constants may assign enormous perturbations to the ele¬ 
ments while at the same time the actual perturbations may be very 
insignificant. 

The original researches of M. Encke on this subject are con¬ 
tained in Nos. 791,792, and 814 of the Astronomische Nachrickten . 
It may be mentioned that a similar method had been already pro¬ 
posed by Mr. G. P. Bond, of Cambridge, U. S. (Memoirs of the 
American Academy of Arts and Sciences , vol. iv. part i. p. 189). 
M. Encke, however, did not become acquainted with this circum¬ 
stance until he had completed his investigation. 


New Method. 
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